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Abstract 

We apply to the XY Z model the technique of construction of integrable models with staggered param- 
eters, presented recently for the XX Z case in jjj. The solution of modified Yang-Baxter equations is found 
and the corresponding integrable zig-zag ladder Hamiltonian is calculated. The result is coinciding with the 
XX Z case in the appropriate limit. 
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1 Introduction 



Recently a technique was proposed which allows to construct a zig-zag ladder type of integrable 
models on the basis of known integrable chain models. It was successfully realised on models with basic 
sl(n) symmetries; sl(2) case (XXZ model) in Jy], sl(3) case (anisotropic t — J model) in [|| and general sl(n) 
case in the || . The main element of this construction is the possibility to stagger the parameters of known 
integrable i?-matrix along chain and time directions. Together with alternation of the anisotropy parameter 
A p], [2| the staggered shift of the spectral parameter u by some additional model parameter 9 was also 
introduced, which caused the appearance of the next to nearest neighbour interaction terms in the resulting 
local Hamiltonian. 

The motivation for considering this type of integrable models is based on the observation that the 
phenomenological Chalker-Coddington model [Q] for edge excitations in the Hall effect after its reformulation 
as a lattice field theory || and calculation of disorder over random phases becomes a Hubbard type model 
with staggered disposition of i?-matrices ||. Therefore it is meaningful to start the investigation of staggered 
models from simple cases. 

In this letter we generalise the construction of the staggered XXZ model and construct a staggered 
integrable XY Z model. The corresponding generalised Yang-Baxter Equations (YBE), which are the 
condition of integrability j?], ||, have a solution. Therefore, one can define a model on a two rung ladder, 
the Hamiltonian of which is presented here. 

2 Basic Definitions and Staggered YBE's 

We start by writing down the modification of basic ingredients of the Algebraic Bethe Ansatz (ABA) 
technique ||] appropriate for our purposes. 

Let us consider Z2 graded quantum Vj tP (v) (with j = 1, N as a chain index) and auxiliary V a>a (u) 

spaces, where p, a = 0, 1 are the grading indices. Consider .R-matrices, which act on the direct product of 
spaces V ata (u) and Vj tP (v), (a,p = 0,1), mapping them on the intertwined direct product of V a ^(u) and 
Vj,p(v) with the complementary a = (1 — a), p = (1 — p) indices 

Raj,ap (U, V) : V a>a ( U ) <g> Vjjv) -> V j>p (v) ® V a ,a(u). (2.1) 

It is convenient to introduce two transmutation operations l\ and L2 with the property i\ = l\ = id for 
the quantum and auxiliary spaces correspondingly, and to define the operators Raj,ap as follows 

Raj,00 = Raj, Raj, 01 = Raji 

Raj,W = R a 2 j, R a j,ll = Ra- 2 - (2.2) 

The introduction of the Z2 grading of quantum spaces in time direction means that we have now two 
monodromy operators T p , p = 0, 1, which act on the space V p (u) = IljLi ^j',p( n ) by mapping it on V p (u) = 

T p (v,u) :V p (u)^V p (u), P = 0,1. (2.3) 

It is clear now, that the monodromy operator of the model, which is defined by translational invariance 
by two steps in the time direction and which determines the partition function, is the product of two 
monodromy operators 

T(v,u)=T (v,u)T 1 (v,u). (2.4) 

The Z2 grading of auxiliary spaces along the chain direction means that the Tq(u,v) and Ti(u,v) mon- 
odromy matrices are defined according to the following staggered product of the R a j(v,u) and R a 2 Av,u) 
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matrices: 



N 



Tl{v,u) = Y\Ra,2j-l{v,u)R L l 2j {v,u) 



N 



T (v,u) = ]jK) 2j „ 1 (v,u)R*%(v,u), 
i=i 



(2.5) 



where the notation R denotes a different parametrisation of the R(v, u)-matrix via spectral parameters v 
and u and can be considered as an operation over R with property R = R. For the integrable models where 
the intertwiner matrix R(v — u) simply depends on the difference of the spectral parameters v and u this 
operation is a shift of its argument u: 



R{u) = R(u), u 



(2.6) 



where 6 is an additional model parameter. We will consider this case in this paper. 

As it is well known in Bethe Ansatz technique M, B, the sufficient condition for the commutativity of 
transfer matrices r(u) = TrT{u) with different spectral parameters is the YBE. For our case we have a set 
of two equations jD 



and 



R 12 (u,v)R^ 3 (u)R 2 3(v) = Ri^(v)Ri3(u)Ri2(u,v) 



R 12 {u,v)R^ 2 (u)R% i {v) = R L £ 2 {v)R L ^{u)R 12 {u,v) . 



(2.7) 



(2i 



3 Staggered XYZ Heisenberg chain 

The ii-matrix of the ordinary XYZ Heisenberg model can be represented as 



R(u) 



( a(u) d{u) \ 

b(u) c(u) 

c(u) b(u) 

\ d(u) a(u) ) 



(3.1) 



Inputting this expression of R(u) into the YBE ( |2.7| ) we obtain the following three sets of equations on 
a{u) , b(u) , c(u) , d(u) . The first set is: 



a(u, v)a Ll (u)a(v) + d(u, v)c L1 (u)d(v) = a Ll (v)a(u)a(u, v) + cf 1 (v)c(u)d(u, v), 
a(u,v)b Ll (u)b(v) + d(u, v)d Ll {u)c{v) = b Ll (v)b(u)d(u,v) + c Ll (v)d(u)d(u,v), 
b(u, v)b ix (u)a(v) + c(u, v)d Ll (u)d(v) = a' 1 (v)b(u)b(u, v) + d Ll (v)d(u)c(u, v), 
b(u,v)a Ll (u)b(v) + c(u, v)^ 1 (u)c(v) = b Ll (v)a(u)b(u,v) + c Ll (v)c(u)c(u, v), 



(3.2) 



which are satisfied automatically in case of trivial i\ operation, as it happened in ordinary XY Z model. 
The second set 

a(u, v)b Ll (u)c(v) + d(u, v)d Ll (u)b(v) = c Ll (v)a(u)b(u, v) + b Ll (v)c(u)c(u, v), 
b(u, v)a L1 (u)c(v) + c(u, v)c L1 (u)b(v) = c Ll (v)b(u)a(u, v) + b Ll (v)d(u)d(u, v), 
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a(u,v)c Ll (u)a(v) + d(u,v)a Ll (u)d(v) = c Ll (v)a(u)c(u,v) + b Ll (v)c(u)b(u,v), 

c(u, v)a L1 (u)c(v) + b(u, v)c L1 (u)b(v) = a Ll (v)c(u)a(u, v) + cf 1 (v)a(u)d(u, v), (3.3) 

b(u, v)c L1 (u)c(v) + c(u, v)a Ll (u)b(v) = a n (v)b(u)c(u, v) + cf 1 (v)d(u)b(u, v), 

b(u, v)d bl (u)d(v) + c(u, v )b Ll (u)a(v) = b Ll (v)a(u)c(u, v) + c' 1 (v)c(u)b(u, v), 

is reducing to staggered XX Z models equations Q in case of d(u) = 0. Finally we have a third set of 
equations 

a(u, v)a L1 (u)d(v) + d(u, v)c L1 {u)a(v) = a' 1 (v)d(u)c(u, v) + d L1 (v)b(u)b(u, v), 

b(u,v)b bl (u)d(v) + c(u, v)d L1 (u)a(v) = a bl (v)c(u)d(u,v) + cf 1 (v)a(u)d(u, v), 

b(u, v)d Ll (u)a(v) + c(u, v)b Ll (u)d(v) = b Ll (v)d(u)d(u, v) + c Ll (v)b(u)d(u, v), 

a(u, v)d bl (u)b(v) + d(u, v)b Ll {u)c{v) = a" 1 (v)d(u)b(u, v) + cf 1 (v)b(u)c(u, v), (3.4) 

a{u,v)c il (u)d{v) + d(u,v)a L1 (u)a(v) = 6' 1 (v)b(u)d(u,v) + c L1 (v)d(u)a(u, v), 

a(u,v)d Ll (u)c(v) + d(u,v)b Ll (u)b(v) = a' 1 (v)a(u)d(u,v) + cf 1 (w)c(-u)a(u, i>), 

which represents the non-trivial part of XYZ model. In the XX Z case, when d(u) = 0, this set of equations 
disappears. 

Our aim is now to find a solution of equations ( |3.2[) , ( |3.3[) , ( |3.4[) with non-trivial l\ operation. It can be 
seen easily that by defining 

a Ll (u) = a(u), d(u, v) = a(u, v ), 

b L1 (u) = -b(u), b(u,v) = -b(u,v) (3.5) 

c L1 (u) = c(u), c(u,v) = c(u,v) 
d Ll (u) = —d(u), d(u,v) = —d(u,v) 



the eqs. ( |3.2| ) obeyed identically, while the eqs. ( |3.3| ) and ( |3.4| ) are reducing to the corresponding equations 
of the staggered case. Therefore the ordinary choice of parametrisation of a(it), b(u), c(u) and d(u) via the 
Jacobi elliptic functions 

a{u) = sn(u + rf) = a*" 1 (u), a(u, v) = sn(u — v + rj) = a(u, v ), 

6(u) = snu = —b Ll {u), 6(it, i;) = sn(-u — v) = —b(u, v), 

c(u) = snr] = c Ll (u), c(u, v) = suit] = c(u, v), (3-6) 

d{u) = /csnr/snnsn(n + rj) = —d Ll (u), 
ci(u, v) = ksnr]sn(u — v)sn(u — v + rj) = —d(u, v), 



fulfils the remaining equations fl3.3|) and 

Let us analyse now the second set of YBE's ( |2.8| ). First we conclude immediately from the solution 
( |3.5[) of the previous set of YBE's that the operation ~ is simply coinciding with the operation t\. Then it 
is easy to see, that if we define 12 operation in fl2.8Q as 



R^(u) = R'^-u) , (3.7) 
then the second set (|2.8|) of YBE's coincides with the first set (|2.7|) after additional action on it by i\. This 



means that the relation ( p.7| ) is ensuring the fulfilment of the second set (2.8) of YBE J s. Therefore we have 
a L2 (u) = a(-u), b L2 {u) = -b(-u), c L2 (u) = c(-u), d L2 (u) = -d(-u). (3.8) 
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It is now necessary to emphasise that, as in ordinary case, the parameters A and k defined as 

a 2 + b 2 -c 2 -d 2 , 2 cd , n , 

= 2a~b ' ^ 77 = (3 ' 9) 

are constants (as well as A Ll and k Ll ). They are defining the anisotropy of the model in the z and y 
directions. As one can see, now, due to solution (p.6j), the anisotropy parameter A is staggered along the 
chain and time directions. 

Hence, we have found the solution of staggered YBE's ( |2.7D , Q2.8| ) and can calculate now the Hamiltonian 
of corresponding integrable model. 

4 The Transfer matrix and the Hamiltonian 

Having the solution of graded YBE's one can start the calculation of the monodromy matrix and the 
Hamiltonian. According to formula fl2.4|) the monodromy matrix of the model is 

t£; e) = t c °; e) t x \ £ : f 2 >, e) 



/ p^l «2iV-2*2JV-l / — \ p &2N-2.?2iV-l ( \ ptli-2 f»2]V-l*2JV / A t/2 &2JV-1 32N ( -\\ (A n\ 

K n Q,2N-l Oajv-iiaAr-l W a 0',2N-l b 2N -ik 2N > tt 0,2N c j 2N \ u ) a 0',2N d k 2N \ U > ) ■ 

Lets write now an explicit formula for the i?-matrix in terms of Pauli matrices 
Ro,r% («) = a(u) (a > r ) + (Jq X j) + 6(u) (a X ] + ctq > r j) + 

<u) (a+ la- ] + *Z X + <) + d(u) (<r+ C V+ * + ^ C V" *) , (4.3) 



where and r refer to auxiliary and quantum spaces respectively and 

, = i (a* + o*) , o = \{o±-<j*), a + = i (a- + to») , = ± (a* _ jo») . 

Then transfer matrix is defined as trace of monodromy matrix ( [4,l|) over auxiliary spaces and at the zero 
value of the spectral parameter u is equal to 

jil—i2n (r\ a\ — rpab,i\...i 2N 



A/ 

= tr tr , n fa* r"t;:U )<rt r "Ho)K: + x(o)^^t 2 ; +1 iL(-^)) ( 4 - 5 ) 



r=l 

N 

= tr tr , J| (sn 2 r/(sn 2 r/ - sn 2 0)Po,2r-Po',2r-i) (4.6) 

r=l 

i.e. up to overall multiplier it is becoming an operator that shifts along the chain by two units, i.e. a 
translation operator. Here we have used 

R Q *(()) = snr? • <S*« = snr? • P , r 
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where Po,r is an operator permuting auxiliary (0-th) and quantum (r-th) spaces. We also used the relation 
(unitarity property of R) 

which holds due to following identities for elliptic functions: 

a{u)a{—u) + d(u)d{—u) = sn 2 // — sn 2 u = b(u)b(—u) + c{u)c(—u) 
a(u)d(—u) + d(u)a(—u) = = b(u)c(—u) + c(u)b(—u). 



The functions a(u), b(u), c(u) and d(u) are given by (|3,6| ) Now we can turn to calculation of the 
Hamiltonian 



— logT 

du 



= sn " sn 2 n) 7i™ Pj^Lw ( 4 ' 8 ) 

u=0 



One can see, that the Hamiltonian of the model is reduces to the sum of contributions from the derivatives 
of quartic products of neighbour i?-matrices inside of brackets in ( f4.5[) after the permutation of indices. It 
appears that the differentiation of two terms in the bracket again contributes to the identity operator and 
provides a constant contribution to Hamiltonian. Therefore the non-trivial contribution looks like 

# l %.Wtf££(0)i^^ (4-9) 

After some calculations by use of properties of elliptic functions and ( |3.9| ) one can find the following 
expression for the Hamiltonian (without the constant part) 

N 

T~L = ^J%2r-l,2r,2r+l,2r+2j (4-10) 
r=l 

Sn 2 77(sn 2 7/ - Sn 2 u)7i 2 r-l,2r,2r+l,2r+2 = (4.11) 
- (a(u)a(-u) - c(u)c(-u)) a'(0) [fffr-i^fr+l ~ <4r a Z r +2] 

+ ( a ( u ) ~ [b'(0) - kc 2 (u)d'(0)) [<r£_ 1 o£. +1 + ^2r-l a tr+l ~ °"2r°2~r+2 ~ fT 2r cr 2r+2] 

( a <X> - (c?(0) - /cc 2 (u)fe'(0)) [o- 2 f r _ 1 cr+ +1 + <J^._ 1 c7^. +1 - (r£ r a£ r+2 - (T 2r a 2r+2 ] 

-b(u)c(u)a'{0) [02r{<Ttr+l a 2r+2 ~ a 2r+l< T 2r+2) + ( CT 2r-l CT 2~r " CT 2~r-l<4r) CT 2r+l] 

~^Y~ ( a ( u ) + a (-«)) - kc 2 (u)d'(0)) [(T2r+l(^tr (7 2r+2 ~ a 2r a 2r+2) + ( a 2r-l a 2r+l ~ a 2r~l a 2r+l) a 2r) 

+ ^y- K«) " «(-«)) (b'(0) " kb 2 (u)d'(0)) l<J2r-l{°tr°2r+l ~ °2r°tr+l) + Kr^r+l " ^^+1)^+2] 

+o(«)d(«)a , (0) [^| r (o-J r+1 o-+ +2 - (T 2r+1 a 2r+2 ) - ((J+_ 1 a+ - o 2r _ x a 2r )a% r+l ] 

+ ^TT (<»(«) + - &C 2 (u)6'(0)) [^+1(0-^^+2 ~ a 2r a 2r+2) ~ ( a 2r-l a 2r+l ~ a 2r-l a 2r+l) a 2r) 



2 

+ C M (a(u ) _ a( _ u)) ( d '( ) _ ^ 2 (n)6'(0)) [^(a+a, 



2r+l ^r^r+l) (^r^r+l a 2r a 2r+l) a 2r+2] 

This expression for the Hamiltonian due to next to nearest neighbour interactions can be easier understood 
as one, written on a zig-zag chain. 
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2j + l 



2j + 3 




2j 



2j + 2 



2j + 4 



Figure 1: Zig-zag ladder chain 

Let us introduce two chains consisting of the even and odd sites of the original chain and label them by 
s = and 1 correspondingly. Make now zig-zag rungs as it is shown in the Fig.l and introduce the following 
labelling of Pauli matrices 



(7 



3,s 



a 2j+s, 



0,1 



(4.12) 



Then, substituting the expression ( j3~lj| ) for a(u), b(u), c(u), d(u) to ( 4. 10|) and after some interesting can- 
cellations of some terms due to identities on elliptic functions, one can write the following zig-zag ladder 
Hamiltonian 

(1 - A; 2 sn 2 7 ? sn 2 6»)(sn 2 7 ? - sn 2 6)n jtS = 



(-1) 



-sn 2 6>cm7dm?(l - k 2 sn 4 r]) [oj, a , } +1)a + 0f jS of +ljS - <rj >s (rf +liS 

i ~abc a b c , ~abc a b e n 1 

t e s a j,s a j,s+l (7 j+l,s "i T s a 3,8° 3,8+1" i+l,8i S — U, 1 

where the anisotropic antisymmetric tensors e" bc and rf 6c are defined as follows: 



(4.13) 



and 



:+3- 



.*3-+ 
e 

2-3+ 



3 = e +3 = -e 3+ = e 3 + = if 3 = -e x +3 = -sntfsnrycnr/dnr^l - fc 2 sn 2 r?sn 2 0), 
" = —£7 3+ = sm?sn6>cn6>dn6>(l - sn 4 r/) (4-14) 



-T, 







T, 







-T, 







-3- _ -+3+ _ --3- 



'0 



T 



-f 3++ = f 3 = — f ] f+3 = fj 3 = — A:snr/sn0cnr/dnr/(sn 2 ?/ — sn 2 #), 

(4.15) 



This tensors define the interaction terms of topological character. The XX Z limit (A;=0) of the Hamiltonian 
( [4.1 3D coincides with the corresponding expression found in W\. 
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